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Non-perturbative quantum geometric effects in Loop Quantum Cosmology predict a modifi- 
cation to the Friedmann equation at high energies. The quadratic term is negative definite and 
can lead to generic bounces when the matter energy density becomes equal to a critical value of 
the order of the Planck density. The non-singular bounce is achieved for arbitrary matter with- 
out violation of positive energy conditions. By performing a qualitative analysis we explore the 
nature of the bounce for inflationary and Cyclic model potentials. For the former we show that 
inflationary trajectories are attractors of the dynamics after the bounce implying that inflation can 
be harmoniously embedded in LQC. For the latter difficulties associated with singularities in cyclic 
models can be overcome. We show that non-singular cyclic models can be constructed with a small 
variation in the original Cyclic model potential by making it slightly positive in the regime where 
scalar field is negative. 

PACS numbers: 98.80.Cq,04.60.Pp 
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I. INTRODUCTION 



One of the cornerstone problems in cosmology is that 
of the initial conditions of our universe. In standard big 
bang cosmology the weight of this problem is shifted to 
quantum gravity which is expected to generate ideal con- 
ditions for the genesis of our universe. Generally, it is 
envisioned that our universe started in a hot expanding 
phase from the big bang and entered a stage of inflation 
before becoming radiation and matter dominated. How- 
ever the idea that our expanding universe was preceded 
by a contracting phase of a classical universe has had 
many avatars. This idea, though promising to eliminate 
some of the fundamental problems related with initial 
conditions in standard cosmology and providing an al- 
ternative to inflation, has remained unsuccessful due to 
absence of generic non-singular evolution from the con- 
tracting to the expanding branch. For a flat fc = cos- 
mological model, which is of interest of the present work, 
based on classical general relativity (GR) such an evolu- 
tion is forbidden by powerful singularity theorems, un- 
less one assumes an exotic form of matter which violates 
positive energy conditions. Though quantum corrections 
have been introduced to achieve a non-singular evolution 
between contracting and expanding phases in special con- 
ditions, it has been realized that a generic non-singular 
evolution is difficult to achieve unless one incorporates 
non-perturbative effects of quantum gravity. In absence 
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of any non-perturbative quantum gravitational modifi- 
cations to dynamics, innovative ideas like pre big bang 
string cosmology and Ekpyrotic/Cyclic models 0j S] 
have so far had limited viability 0. 

Loop quantum gravity (LQG) is a leading non- 
perturbative background independent approach to quan- 
tize gravity 5J. The underlying geometry in LQG is 
discrete and the continuum spacetime is obtained from 
quantum geometry in a large eigenvalue limit. The ap- 
plication of LQG techniques to homogeneous spacetimes 
results in loop quantum cosmology (LQC) ^ which has 
led to important insights on the resolution of singulari- 
ties in various situations B • Using extensive an- 
alytical and numerical methods the analysis of the evo- 
lution of semi-classical states for the flat fc = model 
with a massless scalar field has shown that a contract- 
ing semi-classical universe passes through the quantum 
regime bouncing to an expanding semi-classical universe 
[Tol [ill fl^ . Thus within the context of the model con- 
sidered the idea of the non-singular bounce is realized in 
a natural fashion. 

The underlying dynamics in LQC is governed by a dis- 
crete quantum difference equation in quantum geometry. 
However, using semi-classical states we can construct an 
effective Hamiltonian description on a continuum space- 
time which has been shown to very well approximate the 
quantum dynamics 0,0,0]. In particular, we can ob- 
tain effective equations for the modified Friedmann dy- 
namics which can be used to investigate the role of non- 
perturbative quantum corrections. An important feature 
of the dynamics is that the classical Friedmann equation 
is modified to include a term which is relevant in the 
high energy regime. The modified term in the Friedmann 
equation is negative definite implying a bounce (or more 
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generally a turn-around) when the energy density reaches 
a critical value on the order of the Planck density in ac- 
cordance with the results from the quantum evolution. 

An important question which gains immediate atten- 
tion is the viability of the bounce picture for more gen- 
eral matter sources and compatibility of inflationary and 
Cyclic model ideas with LQC. In this regard various ques- 
tions arise: Does the bounce occur for various inflation- 
ary potentials? Do inflationary solutions form the at- 
tractors of the dynamics after the bounce? Since the 
modifications which lead to a non-singular bounce are 
obtained from non-perturbative quantum gravity effects 
is it possible that these can lead to construction of non- 
singular cyclic models using Cyclic/Bicyclic potentials? 
What insights are obtained to alleviate the problems of 
CycHc models? 

This work seeks to answer these questions using the 
effective theory in LQC. The plan of this paper is as fol- 
lows. In the next section we briefly review the effective 
theory of LQC and classify the differences between quan- 
tum bounce generated by loop quantum dynamics and 
the classical bounce originating from general relativity. 
In Sec. Ill we perform a qualitative analysis using phase 
space diagrams to study the dynamics for inflationary 
potentials {(fP' and (j)^) and show that inflationary trajec- 
tories are attractors of the dynamics after the bounce. 
In Sec. IV we apply the qualitative analysis to nega- 
tive potentials of the Cyclic/Bicyclic models. We show 
that the Bicyclic potential can successfully lead to 
non-singular cycles, whereas the original Cyclic potential 
can exhibit cycling of the scale factor but not the scalar 
field. A truly cyclic model can be easily constructed if the 
Cyclic potential becomes slightly positive in the regime 
where the scalar field is negative, as the Bicyclic poten- 
tial does. The non-singular cyclic behavior obtained for 
negative potentials is unavoidable and occurs for generic 
values of parameters. We conclude the paper with a dis- 
cussion and open issues in Sec V. 



II. EFFECTIVE DYNAMICS IN LOOP 
QUANTUM COSMOLOGY 

LQC is a canonical quantization of homogeneous 
spacetimes based on techniques used in LQC. As in LQC, 
the classical phase space for the gravitational sector is de- 
noted by two conjugate variables the connection and the 
triad, which encode curvature and spatial geometry re- 
spectively. Using symmetries of the homogeneous and 
isotropic spacetime the dynamical part of the connection 
is determined by a single quantity labeled c and likewise 
the triad is determined by a parameter p. For the flat 

= model (which we will only consider in this work) 
on classical solutions of general relativity the relations 
between the new variables and the usual metric variables 
is 

c = 7 d , p = (1) 



where 7 is the Immirzi parameter whose value 7 « 0.2375 
is typically constrained by black hole entropy considera- 
tions 0. 

In the Hamiltonian formulation for homogeneous and 
isotropic space-times the dynamical equations of motion 
can be determined completely by the Hamiltonian con- 
straint. Under quantization, the Hamiltonian constraint 
gets promoted to an operator and the quantum wave- 
functions are annihilated by that operator. It is thus at 
the level of the Hamiltonian constraint that modifications 
due to LQC will appear and from the modified Hamilto- 
nian constraint, the effective Friedmann constraint will 
be derived. Classically in terms of the connection-triads 
variables the Hamiltonian constraint is given by 

Hc1 = -^c2+Hm (2) 

K7^ 

with K = SttG (where G is Newton's gravitational con- 
stant) and TLm being the matter Hamiltonian. The com- 
plete equations of motion are derived from Hamilton's 
equations x = {x^T-Lcx} for any phase space variable x, 
and by enforcing that Tici should vanish. The variables 
c and p are canonically conjugate with Poisson bracket 
{c,p} ~ 7k/3 the use of which gives us relation c = 
from Hamilton's equation for p together with p — a^. 
Substituting these relations into the Hamiltonian con- 
straint ^ and enforcing that the constraint vanishes 
then gives the classical Friedmann equation. 

The elementary variables used for quantization in LQC 
are the triads and holonomies of the connection: hi (/i) = 
exp(^cri) = cos(/xc/2) -I- 2 sin(/ic/2) t^. Here is re- 
lated to Pauli spin matrices as = —iai/2 and /i is 
the eigenvalue of the triad operator p. The gravitational 
constraint is expressed in terms of these elementary vari- 
ables and quantized. On quantization the Hamiltonian 
constraint leads to a discrete quantum difference equa- 
tion whose solutions are non-singular 0, Q • 

The analysis of the quantum Hamiltonian using semi- 
classical states belonging to the physical Hilbcrt space re- 
veals that on backward evolution of our expanding phase 
of the universe, the universe bounces at a critical den- 
sity into a contracting branch. By finding expectation 
values of the Dirac observables we can investigate the 
differences between quantum and classical dynamics. It 
turns out that quantum geometric effects become dom- 
inant only when the energy density (p) of the universe 
becomes of the order of the critical density pcrit EH 0| 
and classical general relativity is a good approximation 
to the quantum dynamics when p <C Pmt ■ 

The quantum dynamics of LQC can be well approxi- 
mated by an effective description. Several propos als have 
been made to derive this effective theory [Ij, [13, [H, [13 
and the result is that the dynamics can be described 
in terms of an effective Hamiltonian constraint which to 
leading order is \2,Q{ 

3 

Heff = a sin^(/[ic) -I- Hm ■ (3) 
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where /i = y 3\/3/2|/i| 1Z| and /i is the eigenvalue of the 
triad operator p in the quantum theory Q : 



pi 
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(4) 



Here Ipi = \JhG denotes the Planck length (in units 
of speed of light equal to unity). In this work we will 
be mainly interested in the matter Hamiltonians corre- 
sponding to a massive scalar field with momentum 110 
and potential Vi^i^Y 



(5) 



The energy density and pressure for the scalar field are 
given by the same expressions as those classically, i.e. 



\i^^-V{<^) (6) 



where we have used Hamilton's equations to get 
H^/p'^/^. It is then straightforward to find the Klein- 
Gordon equation using Hamilton's equations for </> and 
H^, which turns out to be of the same form as the clas- 
sical equation and therefore the scalar field satisfies the 
stress-energy conservation law: 



+ 3-(p0 +P0) = 



(7) 



The modified Friedmann equation can be found by us- 
ing Hamilton's equations for p 

, , 7K 9 2a _ ,- \ ,c\ 

p = {p,7tofr} = — Tr-^i^cS = — sm(/ic) cos(^c) (8) 
3 ac 7/i 

which on using Eq. iQJ implies that the rate of change of 
the scale factor is given by 



a = sin(/xc) cos(/j,c) 
7/i 



(9) 



Furthermore, the vanishing of the Hamiltonian constraint 
implies 



sin^(/ic) 



2,-i2 
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(10) 



Squaring equation ^ and plugging in H1U() yields the 
effective Friedmann equation for the Hubble rate H — 
a/ a 



(11) 



3 \ Pcrit 

with the critical density given by 

V3 



Pcrit 



167r^7'^ 



Ppl 



(12) 



where ppi = \/{hG'^) is the Planck density. The non- 
perturbative quantum geometric effects are manifested 
in the form of a modification of the Friedmann equa- 
tion. Since the modified term is negative definite, when 
P — Petit the Hubble parameter vanishes and the uni- 
verse experiences a turn-around in the scale factor. For 
P ^ Pcrit, the modifications to Friedmann dynamics be- 
come negligible and we recover standard Friedmann dy- 
namics. Note that origin of pcrit is purely quantum, since 
as Ppl cx 1/ft, pcrit ^ 00 in the classical limit /i — > 0. 

Interestingly, similar p^ modifications have also been 
obtained in Randall-Sundrum braneworld scenarios 
based on string theory but those come with a positive 
sign in the Friedmann equation and a bounce is not pos- 
sible. If one assumes the existence of a time-like extra 
dimension then the sign can become negative as in LQC 
[2^ 1 (for an analysis of dynamics of the brane- world model 
with two time-like dimensions see Ref. '25^ and "26]). A 
detailed comparison of the effective dynamics of LQC and 
braneworld scenarios shows interesting features which in- 
clude a dual relationship between their scaling solutions 



III. QUANTUM AND CLASSICAL BOUNCES 

Having established the equations for the effective 
Friedmann dynamics we now wish to classify the condi- 
tions under which effective dynamics given by the mod- 
ified Friedmann equation Hll(l implies a turn-around; 
namely either a bounce or re-collapse. As discussed above 
when the energy density of the universe reaches the crit- 
ical value pcrit, the Hubble parameter becomes zero and 
a bounce/re-collapse is expected. Since we will also con- 
sider the possibility of the scalar field potential being 
negative we also discuss the case where the matter en- 
ergy density vanishes which also leads to a turn-around. 

More precisely from Eq. we can determine the 
conditions under which such a turn- around can occur; 
i.e., when d = 0. It is clear from that equation that 
there are two possibilities, namely when pc — (rt -I- ^)7r 
and pc = rriTT for any integers n, m. The energy density 
can be computed for both cases and is equal to pcrit for 
the former and zero for the latter. We thus label turn- 
arounds occurring due quantum geometric effects when 
pc = {n+ i)7r as quantum turn-arounds and those which 
occur when energy density vanishes for pc = mir as clas- 
sical turn-arounds. 

To determine whether or not a bounce or recoUapse 
occurs we need to go further and calculate the second 
time derivative of the triad (which is proportional to d 
at the turn-around), where a negative value indicates a 
recollapse and a positive value a bounce. Taking the time 
derivative of Eq. (ISJ and evaluating when p equals zero 
gives 
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p \p=o — — — cos{Zpc)c . 



(13) 
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To determine the sign of p we therefore need to know the 
time derivative of c which is determined from Hamilton's 
equations c = {c, Tieff} leading to 

■ sin(/ic) cos(/ic) — 



pi 



sm (//cj ^ 



3> 



3 

1 dUm 
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(14) 



where the first term vanishes for both quantum and clas- 
sical turn-arounds. We now need to consider separately 
the quantum and classical turn-arounds and to specify 
the form of matter. We first consider the quantum turn- 
arounds. 



A. Quantum Turn-arounds 

The quantum turn-arounds as stated above are a di- 
rect consequence of the discrete quantum geometry ef- 
fects of LQC. From the fact that /Ic = (n -I- i)7r we have 
sin(/2c) = 1 and cos(/ic) — which implies the matter 
energy density satisfies p = Pcrit- Under these conditions 
the equation (|14|l for the triad acceleration becomes 



P \p=0 = up Pent - ^ 



2 1 dT-L 



M 



3^ dp 



(15) 



Let us first consider the simplest case of matter with 
constant equation of state w with energy density p oc 
^-3(1+J«) ^ p-3/2(l+t«) ^^^^^ implies Hm oc p-3»/2. It is 

then straightforward to calculate p giving 



p \p=0 = Kppcnt{l + w) 



(16) 



and thus we can make the identification that for w < 
— la recoUapse can occur and for w > — 1 a bounce 
can occur. This can be understood in a simple fashion. 
Matter with equation of state w < —1 violates the null 
energy condition and hence has increasing energy density 
as the universe expands. Thus when the energy density 
hits the quantum critical value a re-collapse occurs and 
the universe begins contracting. In the classical universe 
the energy density of such matter can grow to infinite 
values and even lead to a singularity (for example a big 
rip singularity encountered for phantom field which has 
a negative kinetic energy '28^). In the context of LQC 
the implication is that such a singularity does not occur 
as the e nerg y density remains bounded and the big rip is 
avoided |23|- On the other hand matter with w > —1 has 
increasing energy density as the universe contracts and 
therefore a quantum bounce will occur when the critical 
density is reached. 

Let us now turn to the case of a massive scalar field. 
Here we consider general forms of the potential which 
can be negative valued and phantom models where the 



kinetic term is negative definite. The general form of the 
matter Hamiltonian is 

7iM = ±^p-'/'H2 +//2y(^) (17) 

where the ± indicates normal and phantom matter re- 
spectively. With this form wc can calculate which 
gives 



OHm 3 / 1 



T^P^'^'nl+p'/'V{4,)] (18) 



dp 2p\ 2 



and we can simplify this by noting that when p = pcrit we 
have ±ip^3/2jj2 _ p3/'2^_y _|_ p^^.j.) whence p simplifies 

to 



p |p=o = 2kp (pcrit - V) 



(19) 



The implication of this is that if the scalar field potential 
is less that the critical density at the turn-around then 
a bounce occurs, otherwise if the potential is larger a 
recollapse occurs. In this form we can analyze the turn- 
around behavior for the different types of scalar field and 
potential. 

We start with a normal scalar field where the kinetic 
term of the matter energy density is positive definite. Be- 
cause of this when the total energy density is equal to the 
critical value at the turn-around the potential must be 
less than than pcrit- Thus for a normal scalar field, only 
a bounce can occur in a quantum turn-around. This can 
be understood intuitively from the previous paragraph 
given that an ordinary scalar field can not have w < ~\. 
Since we have not made any assumptions on the form 
of the potential even negative potentials will not lead to 
a quantum recollapse (though as we shall see a negative 
potential can lead to a classical recollapse). 

Considering now phantom models with a negative def- 
inite kinetic term we find the opposite situation. First 
a positive potential is required in general from the effec- 
tive Friedmann equation to get classical behavior of the 
scale factor. Furthermore since the kinetic term is nega- 
tive definite when the matter density equals the critical 
density, the potential is larger than the critical density. 
Hence p <0 and the phantom field behaves like a matter 
with w < —1 and a quantum recollapse can occur. This 
implies generically that a phantom field can not have a 
quantum bounce in LQC. 



B. Classical Turn-Arounds 

The classical turn-around, characterized by pc = rmr 
and p = 0, is not a feature limited to LQC but would 
occur even in the classical theory. For standard perfect 
fluids with fixed equation of state this does not occur as 
the energy density is always non-zero. However, if the 
scalar field potential is negative or a phantom field is 
present this condition can be met. At this turn-around 
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cos(/ic) = 1, sin(/ic) = which gives us 



P\p=0 



-up 



1/2 



M 



3 dp 



(20) 



and we can again use p = to relate ^\p '^^^H^ — p^l'^V 
giving 



p |p=o = InpV . 



(21) 



In the case of a normal scalar field the potential must 
be negative in order for the energy density to become zero 
hence p < at turn-around and a recoUapse is possible. 
In the case of a phantom model we require the poten- 
tial must be positive and hence a classical bounce can 
occur. Furthermore, in the course of dynamics it is pos- 
sible that both a quantum and classical turn-around can 
occur. This is precisely what can happen when we con- 
sider negative potentials as the universe cycles between 
classical recollapses and quantum bounces. 



IV. QUALITATIVE ANALYSIS OF BOUNCE 
WITH INFLATIONARY POTENTIALS 

The analysis of the previous section has established an- 
alytically the turn-around behavior if the conditions for a 
quantum or classical turn-around are met. However the 
analysis does not indicate whether the turn-around is a 
feature of the dynamics given generic initial conditions. 
To go further we will use phase portraits to analyze the 
space of solutions to the dynamics. To illustrate the tech- 
niques and as a simple non-trivial model we first consider 
an ordinary scalar field with (j)'^ and cj)'^ potentials con- 
ventionally used for chaotic inflationary scenarios. In the 
previous section we showed that only a quantum bounce 
is possible. We are particularly interested in the question 
if inflation in the post bounce expanding universe is an 
attractor in the space of solutions. 

It is useful to first obtain H equation using Friedmann 
equation Hll|) and the conservation law 0: 



H : 



iP + P4>) 1 



2p 

Pcrit 



(22) 



where p = © . 

Already from 1)11(1 and iQ it is evident that both energy 
density and Hubble parameter are bounded 



< P < Pcrit, 



H{p = Pcrit) =H{p = {))={) 



1 

2 Per it 5 



(23) 
(24) 
(25) 
(26) 



For the scalar field with non-negative potentials, the only 
nontrivial case is a quantum turn-around, i?(/Ocrit) = 



and it always represents a bounce, since p + p^ = (j)'^ > 
and 7? > in (|22|) . Therefore, given any positive definite 
potentials all solutions of the cosmological equations (UJ 
and (|ll|l are non-singular. As an example, let us discuss 
the simplest case of a massless scalar field. The solution 
of O can be written as 



— V^Pcrit ( — 



(27) 



Substituting this solution in Eq. Hll|l gives the analytic 
solution for Hubble parameter 



(28) 



Clearly this solution is devoid of any singularity as 
soon before and after the bounce the Hubble parameter 
reaches its maximum value. 

When the effective potential of the scalar field does not 
vanish analytical solutions are difficult to find. To under- 
stand the qualitative behavior of solutions it is useful to 
build phase portraits. The dynamical phase space is four 
dimensional (consisting of c{t) , p{t) , (l){t) , (j){t)) however 
the vanishing of the constraint can be used to fix one of 
them given the other three. Therefore a complete phase 
portrait would be three-dimensional, however for simplic- 
ity we will display two dimensional portraits consisting of 
(/> and (f). Since the Hubble parameter only changes sign 
at a classical or quantum turn-around each phase portrait 
drawn will consist of either an expanding or contracting 
phase. The quantum and classical turn-arounds will be 
represented as boundaries in the phase portrait (which 
we will represent as solid lines in the figures) and any 
trajectory that reaches the boundary will continue into 
the corresponding phase portrait with opposite sign of 
the Hubble rate thus indicating the turn-around. In this 
sense the expanding and contracting phase portraits are 
to be glued together at the boundaries. 

In order to obtain equations for corresponding dynami- 
cal system we can express the Hubble parameter in terms 
of (j) and (p using the effective Friedmann equation and 
plug into the conservation law Eq. to get 



dv 



V 



1 



ir(* + '')l} 



(29) 

where the sign "-" corresponds to expansion while the 
sign gives contraction. Using this equation we can 
obtain the phase space portraits [s^] for the massive 



scalar field V{(j)) = m (p /2 with 



m 



0.2 and self- 



interacting nonlinear scalar field with V{4)) — Xcj)^ /A with 
A = 10^^. As we have discussed in the previous section, 
in this case only quantum bounce is possible which occurs 
when p = Pcrit, that is 



+ v{cp) = pc 



(30) 



Regions close to the outer boundary correspond to high 
energy density with quantum effects dominant whereas 
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the low energy limit occurs near the origin in the (f), (j) 
plane. 




Figure 1: Phase portrait for massive scalar field 0^ potential. 
Dashed curves represent GR case and solid curves shown LQC 
case. Dots correspond to points where initial conditions for 
numerical solutions were given. The direction of phase tra- 
jectories is from the boundary to the origin. Ifere m = 0.2. 

The results for the two potentials in the expanding 
phase are shown in fig. ^ and fig. [3 respectively. The 
phase space trajectories obtained via LQC dynamics are 
compared with those of GR obtained in Ref. jSlj which 
are given as the dotted lines in the figures. All phase tra- 
jectories are directed towards the origin which is the only 
particular point in the finite region of phase space having 
the character of a stable node. The scalar field experi- 
ences numerous damped oscillation around the origin in 
agreement with classical GR. The agreement between the 
two theories is indicative of the fact that p <C Pcmt near 
the origin. In the region near the outer boundary the 
LQC trajectories deviate radically from the GR case and 
the existence of the boundary is a fundamental difference 
between LQC and GR 32] (see analogous discussion for 
gauge theories of gravity case in [ssLls^ ). 

A key feature of these phase portraits is the inflation- 
ary separatrix (see details for e.g. in Refs. 0,1^ 0|) 
which appears in the figures as the curve in the second 
and fourth quadrants to which the trajectories are at- 
tracted before undergoing oscillatory behavior. It is here 
that inflationary behavior occurs and it was shown in Ref. 
[sif that most solutions for (jp' and 0"* potentials within 
GR tend to this separatrix in the expansion stage. This 
led the authors of Ref. [sj to the conclusion that in- 
flation is generic within these models. This separatrix is 




Figure 2: Phase portrait for the self-interacting scalar field 
cj)'^ potential. Trajectories for LQC are shown by solid curves 
where as dashed curves show GR case. Here A = 10~^. 



also present in the LQC case and is qualitatively the same 
although its position slightly changes with respect to the 
case of GR especially in the vicinity of the boundary. 

However, in LQC the trajectories do not begin at the 
outer boundary but are preceded by a contracting phase 
before the bounce. The phase portraits can be obtained 
for the contracting phase by reflecting the expanding por- 
traits with respect to the vertical axis and inverting the 
direction of the trajectories. Now all cosmological solu- 
tions start in the unstable node at the origin where the 
universe is contracting and the scalar fleld is behaving 
as a driven oscillator (anti-friction like behavior), reach 
the outer boundary giving a bounce, return to the in- 
flationary separatrix in the expanding phase and finally 
approach the stable node at the origin where the scalar 
field experiences damped oscillations. We show the com- 
plete phase portrait of massive scalar field within LQC 
in figOlby superimposing the contracting and expanding 
phases. The dashed curves correspond to the contract- 
ing phase and full curves describe the expanding phase. 
The thick curve gives a concrete example of the closed 
phase trajectory which starts at the origin, comes to the 
boundary, experiences a non-singular bounce and comes 
back to the origin. 

Due to the fact that infiationary separatrix is also 
present in the case of LQC we come to the conclusion 
that infiation is as common feature of cosmological so- 
lutions within LQC as it is within GR. If one defines a 
measure of solutions on the boundary (at the bounce) 
it turns out that only a small fraction of solutions does 
not contain infiationary phase. Since all solutions are 
non-singular, they continue in the past to the contract- 



7 



ing phase. The corresponding separatrix has a repulsive 
nature in contracting phase (i.e. solutions tend to de- 
viate more and more from it in course of cosmological 
contraction), and the probability for solutions to follow 
this seperatix, given equal measure of initial conditions 
at the bounce, is small 35]. In fact this separatrix was 
shown to be exponentially unstable within GR, meaning 
that solution which is initially close to it deviates expo- 
nentially from it in course of time js^. Therefore, if all 
initial conditions have equal probability at the boundary, 
the probability to get non-inflating solution is exponen- 
tially suppressed. An example of the phase trajectory in 
figEI shows that solution can not follow repulsive separa- 
trix on contracting stage but nevertheless has inflationary 
stage during expansion. 



V. QUALITATIVE ANALYSIS OF BOUNCE 
WITH NEGATIVE POTENTIALS 

Negative potentials have been used in Cyclic models to 
construct alternatives to inflationary scenarios. It is en- 
visioned that seeds for the generation of cosmic structure 
originated in the contracting phase of the universe which 
preceded the current expanding phase. A basic difficulty 
in these scenarios is the singular nature of transition be- 
tween two phases which has been extensively studied in 
Ref. As we noted in Sec. II, a classical bounce is 

not possible with negative potentials for a normal scalar 




field, whereas a bounce mediated by loop quantum cor- 
rections is possible. We now analyze the nature of this 
bounce for different negative potentials. 



A. Massive scalar field 

For simplicity we start with massive scalar field with 
a constant negative potential of form 



(31) 



where Vb > 0. Without loss of generality for the nu- 
merical simulations we use Vb = 0.1 and m — 0.2 in this 
subsection. For this potential there exists an inner and an 
outer boundary in phase portraits. The inner boundary 
arises due to classical recollapse, a feature shared with 
classical GR The outer boundary represents the 

quantum bounce and is absent classically. The presence 
of the two boundaries in LQC opens a novel possibility 
to construct cyclic models where the universe undergoes 
a series of expansion and contraction phases. 




H 




Figure 4: Behavior of 4>, (j) and H is shown for the mas- 
sive scalar field potential 1)3 1^ in LQC. The field starts rolling 
down the negative potential from positive value of <j) and turns 
around due to loop quantum effects reaching the value where 
it started from. The scale factor also bounces as is evident 
from the turn-around in the sign of Hubble parameter. Initial 
conditions are cj> — 0.640, (j) — 5.422; parameters are Vo = 0.1, 
m = 0.2. 



Figure 3: Complete phase portrait for massive scalar field 
within LQC showing both contracting (dashed curve) and ex- 
panding phases (full curve). Thick curve shows a concrete 
example of the phase trajectory. 



Numerical solutions for the time evolution given a set 
of initial conditions are shown in fig0|and the phase por- 
trait is shown in fig[Sl Both the Hubble parameter and 
energy density are bounded subject to the constraints 
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given in H23() and H24() . As can be seen from fig0I the 
scalar field starts from positive values and rolls down 
the potential giving rise to a period of inflation with the 
Hubble parameter positive and almost constant. This 
continues until the scalar field enters the region of neg- 
ative potential and the energy density reaches zero at 
which point the classical recoUapse occurs and H — 0. 
In classical GR the field would continue its motion to- 
wards negative (j) with ever increasing (because of the 
anti-friction term in the Klein-Gordon equation) and fi- 
nally ending in a big crunch singularity. However, in LQC 
when (f) reaches a value such that the energy density of 
the field becomes comparable to Pcrit, the trajectory de- 
viates from the classical one and the magnitude of the 
Hubble parameter starts decreasing and quickly becomes 
zero at p = Pcrit- The universe bounces and immediately 
afterward enters a phase of super-inflation, where H > 
(which has been shown to be a generic feature in LQG for 
P > Pcrit/2 |27j). Since the universe expands quickly af- 
ter the bounce, (p decreases and p becomes smaller than 
/3crit/2. The phase of super- inflation ends and Hubble 
rate starts decreasing. The field is now negative valued 
and climbs up the potential only to turn-around, undergo 
slow-roll inflation and repeat the process in the opposite 
direction. The cycling continues indefinitely which is ev- 
ident from the phase portrait in figlH 

As can be seen from fig01 the oscillations of the dy- 
namical variables are asymmetric. This is due to the fact 

4) 



Figure 5: Phase trajectory of the solution shown in figB] The 
inner boundary corresponds to classical bounce (p — 0) and 
the outer boundary to quantum bounce (p — pait)- The 
trajectory formally corresponds to infinite asymmetric oscil- 
lations, although is shown for a finite time interval. 



that during expansion the inflationary separatrix (which 
is still present, like in the cases discussed in Sec. HI) has 
an attractive nature, while during contraction it is repul- 
sive one. Therefore the expanding universe spends more 
time in the inflationary phase as opposed to the time 
spent in the contracting phase and the universe experi- 
ences overall net growth in the scale factor. The oscilla- 
tions of the scalar field together with a monotonic Hubble 
parameter, which are typical for usual chaotic inflation- 
ary potentials, occur only if the inner boundary is close 
enough to the origin (corresponding to a small value of 
Vq) allowing for the scalar field to oscillate before being 
reflected by the inner boundary. 

As is evident from the phase portrait in figEl given 
the parameters chosen these oscillations do not occur. 
The presence of the two boundaries leads to the only 
one global direction of phase trajectories, the clockwise 
direction seen from fig|Sl 

B. Cyclic potential 

Let us now consider the cyclic potential in the form Q 

F = Vb(l-e~'"'^)exp(-e-"'*). (32) 

In the Gyclic model, the value Vq is chosen so as to give 
the correct magnitude of the current cosmic acceleration 
and the other parameters are typically constrained to 
give the correct amplitude for density perturbations. For 
the purposes of qualitative analysis and without any loss 
of generality, we choose the following parameters: a = 
0.3\/8^, cj = 0.09^8^, and = 0.1. The minimum of 
the potential is displaced from the origin and the inner 
boundary is situated to the left of the origin. The phase 
diagrams for expansion and contraction stages with this 
potential are represented in figEl 

Here, as in the previous case, both a classical recol- 
lapse and quantum bounce can occur. Again, there is a 
separatrix in the phase space located to the right from 
the inner boundary. As in the previous cases during con- 
traction phase trajectory is directed towards the outer 
boundary and during expansion it is directed towards 
the inner boundary. 

The evolution of the Hubble parameter in LQG with 
the Gyclic model potential can be understood from fig. 
[3 As in the previous case we investigate the dynamics 
when the field starts from positive </) where the potential 
is positive. The field rolls down towards negative val- 
ues of (j) and upon entering the region of negative poten- 
tial its kinetic energy cancels the potential energy lead- 
ing to -ff = and the classical recollapse. The universe 
starts contracting, heading towards a big crunch but ex- 
periences a quantum bounce when p — pcrit- After the 
bounce, the universe briefly enters the phase of super- 
inflation before the Hubble parameter starts decreasing. 
However, unlike the previous case the potential is still 
negative valued and the field can not turn around and 
continues to cf) = —oo. Another round of collapse and 
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bounce occurs, but this behavior terminates and the uni- 
verse approaches the state H ^ with (j) ^ and ^ 
at i — > oo. Thus in the case of cyclic model potential 
though the scale factor bounces and universe escapes the 
big crunch singularity, the dynamics does not lead to cy- 
cles. For such cycles to exist the scalar field must turn 
around and return to the region of positive (/). To achieve 
this the potential must become positive at some point 




Figure 6: Phase diagrams for contraction (above) and expan- 
sion (below) stages with the cyclic potential H32|l . Red curve 
shows repulsive (attractive) inflationary separatrix. Here 
a = 0.3-\/87r, lo = 0.09\/87r, and Vo = 0.1. Outer boundary 
corresponds to p = pcrit and inner boundary to p = 0. 
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Figure 7: Behavior of 0, <j) and H with the Cyclic model 
potential 13211 . Initial conditions for this solution are: tfi = 
16/3, (f) = 1.066. Parameters are: a = O.SVStt, uj — O.OQv^Svr, 
and Vo = 0.1 



for negative (j) as was the case for the quadratic potential 
with negative constant. Since the Cyclic model potential 
is negative for all values of < and approaches zero 
asymptotically, the possibility of such a turn-around in 
(j) does not exist. 

The structure of the phase space shown at figEl sug- 
gests that all solutions begin with cf) — —oo and (p > 
and end at </) = — oo with (f> < 0. The boundaries define 
the overall direction for the phase space flow which is 
clockwise, like in the case of massive negative potential, 
but the outer boundary is not closed, so the scalar field 
escapes to negative infinity. This is illustrated at figlHl 
with a particular solution. We have set initial conditions 
in a way to have positive (j) at the beginning with the 
purpose to illustrate the dynamics with cyclic potential 
within LQC. However the evolution can be traced con- 
tinuously to the past with the field having large negative 
values. 

The problem of non-singular bounces for C ycli c po- 
tential in LQC has been investigated earlier [331, but 
by considering modifications to only matter part of the 
Hamiltonian. The analysis of Ref. j33| was able to show 
the existence of non-singular bounce in scale factor for 
Cyclic model for some choices of parameters but not in 
general. A limitation of that analysis was the problem 
of energy density becoming super-Planckian as the field 
rolls down the Cyclic potential [s^- As we understand 
from the present analysis, incorporation of modifications 
to the gravitational part in the effective Hamiltonian au- 
tomatically solve this problem, as the energy density is 
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bounded above by a critical value. 



4) 




Figure 8: Phase portrait of solution shown at figEI for wider 
time interval. Here again we draw only finite time evolution 
of the solution. Inner boundary corresponds to p = and 
outer to p = Pcrit- 



C. Bicyclic potential 

The Bicyclic potential was discussed in Refs. |3, QD| 
as a better choice over the Cyclic potential in the at- 
tempt to construct cosmologically viable cyclic models. 
In fact there exists the crucial difference between the 
cyclic and inflationary models. The principle difficulty 
of the former is the cosmic singularity, but in a much 
more drastic way than appears in standard inflation. In 
contrast to inflation, within cyclic models perturbations 
are created at the contraction stage and cannot be traced 
through singularity. Instead, in our case, given regular 
background dynamics perturbations can in principle be 
evolved through the regular bounce once these are incor- 
porated in LQC. At present we are interested in compar- 
ing its phase portraits with the Cyclic potential. 

The form of the Bicyclic potential we will discuss is 

V ^Vo[l- Acosir^{(j))] , (33) 

where A > and Vq > 0. For qualitative analysis we 
choose A = 10 and Vq — 0.03. The phase diagram for 
the expansion stage is represented in figlHl As in the pre- 
vious cases there exist inner and outer boundaries corre- 
sponding to classical recoUapses and quantum bounces. 
As can be seen from the phase portrait the outer bound- 
ary becomes a horizontal line far from the origin given 



by the condition 4>b = ^V^Pcrit since for large values 
of (f) the potential is small. The inner boundary is a 



<i> 




Figure 9: Phase diagram for expansion stage with the Bicyclic 
potential 1331 1. Red curve again shows attractive inflationary 
separatrix. Inner boundary refers to p = and outer to p = 
Pcrit. Here A = W and Vo = 0.03. 

closed curve surrounding the origin crossing the (j) axis 
at points where the potential (|33|l vanishes, namely when 
(j) = ±arccosh(A). 

The phase portraits for the contracting branch can be 
obtained by reflecting the phase diagram in figEI with 
respect to the vertical axis and inverting the directions 
of phase trajectories. Example of phase trajectory for 
bicyclic potential is represented at figllOl 

A difference between the Bicyclic model and the previ- 
ous models is that there is a closed inner boundary but an 
open outer boundary. However as in the case of massive 
scalar field with negative constant there is a preferred 
clockwise direction for phase trajectories. In contrast to 
the Cyclic model, the behavior cycles both in terms of 
the Hubble expansion and the scalar field. Given initial 
conditions between the boundaries, the solution is at- 
tracted to the inflationary separatrix until reaching the 
inner boundary. After recoUapsing the separatrix is re- 
pulsive and the solution is driven to the outer bound- 
ary leading to a quantum bounce and the process starts 
anew. Thus, with the Bicyclic potential there is an infi- 
nite number of inflationary stages within a given solution. 
The time variation of (j), (jj and H is qualitatively similar 
to that for the massive scalar field shown in fig0]with a 
bounce in both the scalar field and the scale factor. The 
Bicyclic model thus leads to a non-singular cyclic model 
in LQC. Again the crucial distinction between the Cyclic 
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<i> 




Figure 10: Phase portrait with the Bicyclic potential H33^ for 
a finite time interval. Outer boundary corresponds to p = pcrit 
and inner boundary refers to p = 0. 

and Bicyclic potential is the fact that the potential is 
positive for both large negative and positive values of (j). 
This cycling in the Bicyclic model is responsible for the 
better viability of the latter. 

VI. DISCUSSION AND CONCLUSIONS 

The non-perturbative quantum geometric effects lead 
to a modification with a negative sign in the Fried- 
mann equation in LQC. Since the correction term is neg- 
ative definite it can lead to a quantum bounce in the 
high energy regime when loop quantum modifications are 
dominant. We have investigated the qualitative details 
of this bounce for inflationary and negative potentials. 
The existence of the outer boundary given by ^ in the 
phase portraits guarantees non-singular behavior of all 
the solutions with a scalar field in LQC with any kind 
of potential. For negative potentials the inner bound- 
ary also appears corresponding to the classical recoUapse. 
The presence of the two boundaries for negative poten- 
tials leads to the possibility of solutions having cyclic 
behavior. We have also briefly reported on the nature of 
quantum turn-arounds for a phantom field. 

The massless scalar field gives a good example of one 



feature of cosmological solutions within LQC, namely the 
non-singular bounce. In effect, at least solutions with 
(/)^ and 0^ potentials can be well approximated near the 
bounce by the corresponding exact solutions for a mass- 
less scalar field, because most of them approach the outer 
boundary of phase space with kinetic energy much larger 
than the potential energy, thus the potential can be ne- 
glected. For these potentials we showed that inflationary 
trajectories are attractors of dynamics after the bounce. 
In the same way a massive scalar field with negative po- 
tential represents a good example of cycling between ex- 
pansion and contraction accompanied with inflation, ex- 
hibiting the main features of the more complicated po- 
tentials. Our analysis shows that there is little quali- 
tative difference between the massive scalar field with 
negative constant and Bicyclic potentials, although their 
functional form is very different. In both cases cyclic 
solutions can be found with an infinite number of infla- 
tionary stages. 

Though the problem of the big crunch can be over- 
come with potential of the Cyclic model, a limitation 
remains in the lack of a turn-around of (j) from the neg- 
ative side of the potential. In the Cyclic model this in- 
stant is envisioned as the collision of two branes. If the 
effective potential between branes can become positive 
prior to collision (arising from higher order perturba- 
tive string effects/or non-perturbative corrections), then 
a turn-around of the field can occur and a viable cyclic 
model can be constructed. On the other hand, the Bi- 
cyclic potential provides the possibility for the field to 
turn around and may be used for further development of 
cyclic models. 

The qualitative analysis reported in this work 
has shown that both inflationary paradigm and the 
Cyclic/Bi-cyclic scenarios may be incorporated in LQC. 
Our analysis has not touched on the aspects which can 
quantify the viable parameter space of these models with 
respect to observational data. Further, we have not ex- 
plored the way LQC may shed insights on the problem of 
evolution of perturbations through the bounce. Work on 
including perturbations in the quantum theory and de- 
riving effective equations is in progress, which will even- 
tually address the viability of these ideas with respect to 
observations. 
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